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We apply a recently constructed model of analytic QCD in the Operator Product Expansion 
(OPE) analysis of the r lepton decay data in the V+A channel. The model has the running 
coupling Ai(Q 2 ) with no unphysical singularities, i.e., it is analytic. It differs from the corresponding 
perturbative QCD coupling a(Q 2 ) at high squared momenta \Q 2 \ by terms oc (1/Q 2 ) 5 , hence it does 
not contradict the ITEP OPE philosophy and can be consistently applied with OPE up to terms 
of dimension D = 8. In evaluations for the Adler function we use a Pade-related renormalization- 
scale-independent resummation, applicable in any analytic QCD model. Applying the Borel sum 
rules in the Q 2 plane along rays of the complex Borel scale and comparing with ALEPH data of 
1998, we obtain the gluon condensate value ((a 3 /ir)G 2 } = 0.0055 ± 0.0047 GeV 4 . Consideration 
of the D = 6 term gives us the result {0 { ^ +A) ) = (-0.5 ± 1.1) x 10~ 3 GeV 6 , not incompatible 
with nonnegative values. The real Borel transform gives us then, for the central values of the two 
condensates, a good agreement with the experimental results in the entire considered interval of 
the Borel scales M 2 . In perturbative QCD in MS scheme we deduce similar result for the gluon 
condensate, 0.0059 ± 0.0049 GeV 4 , but the value of D = 6 condensate is negative, (Og V+A ') = 
(—1.8 ± 0.9) x 10 -3 GeV 6 , and the resulting real Borel transform for the central values is close to 
the lower bound of the experimental band. 

PACS numbers: 11. 10. Hi, 11.55.Hx, 12.38.Cy, 12.38.Aw 



I. INTRODUCTION 



Perturbative QCD (pQCD) in the usual renormalization schemes, such as MS, has the peculiar property of resulting 
in a running coupling a(Q 2 ) (= a s {Q 2 ) /it) which has singularities in the complex Q 2 plane outside of the negative 
semiaxis, often at positive Q 2 = —q 2 > 0. This implies that spacelike physical quantities T>(Q 2 ) evaluated in terms 
of a{Q 2 ) have the same type of singularities, contravening the general principles of (causal and local) quantum field 
theories 0, HJ . These problematic aspects of pQCD were addressed in the seminal works of Shirkov, Solovtsov et al. 0- 
[E[ who constructed, as an alternative, a QCD model which has no such unphysical singularities. Their idea was to 
keep the discontinuity function p^f^^a) = Im a(— a — ie) unchanged for Q 2 < (c > 0), but eliminating the offending 
singularities at Q 2 > by not including them in the dispersion relation for their coupling. The same principle was 
applied to any other integer power of a(Q 2 ). This approach was called Analytic Perturbation Theory (APT); in the 
present context we call it Minimal Analytic (MA) QCD, due to the unchanged ( per turbative) singularities along the 
Q 2 < semiaxis. For various applications of MA and related models, see Refs. 043 - 1 However, MA has two aspects 
which, under specific circumstances, may be regarded as inconvenient: 

1. It subestimates (HQ the semihadronic r decay ratio 2 rv, giving r T (MA) « 0.13-0.14 in the strangeless V+A 
channel while the experimental value is 0.203 ± 0.004 (181 Il9j. 

2. At momenta Q 2 > A 2 (where A 2 ~ 0.1 GeV 2 ), the MA coupling „4^ MA) (Q 2 ) differs from the underlying QCD 
coupling a(Q 2 ) by terms ~ (A 2 /Q 2 ). The latter property implies that, in MA, the leading-twist contributions to 



1 In comparison with v2: rewritten and extended Section II; Section III refers to a new Appendix A; the results in Section IV are 
obtained by the standard minimization involving \ 2 i Conclusions include a comparison with other results in the literature for the gluon 
condensate; new references [47] and [58]. To appear in PRD. 

1 For a somewhat related but different approach, which performs minimal analytization of dlna(Q 2 )/dlnQ 2 (and not: a(Q 2 )), see 
Refs. [l5l ]. and for extensions thereof Refs. [id ]. 

2 r r represents the QCD part of the strangeless V+A decay ratio R T (AS = 0); it has the (small) quark mass effects subtracted and is 
normalized in the canonical way: r T = a + 0(a 2 ). In Ref. % the correct value of r T is reproduced in MA at a price of (modifying MA 
by) introducing strong mass threshold effects in the low-momentum regime. 
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physical quantities contain power terms ~ (A 2 /Q 2 ) which are of ultraviolet origin. This is not in accordance with 
the ITEP interpretation of OPE [13, which states that the higher-dimension (higher-twist) terms ~ (A 2 /Q 2 ) n 
in OPE are of infrared origin. 



These problems have been addressed in a series of papers [22H24j . In Ref. [22| the following question was investigated: 
Is there an analytic QCD which is simultaneously fully perturbative and reproduces the correct value of r r ? The 
results of Ref. 22] indicate that there is probably no such framework, unless we introduce renormalization schemes 
which make the perturbation series highly divergent starting at terms ~ a 5 . Therefore, in Refs. [2^, [24| a more 
modest goal was pursued: construction of analytic QCD which is not fully perturbative, but addresses the previously 
mentioned two points nonetheless: the correct value of r T is reproduced, and the deviation of the analytic coupling 
Ai(Q 2 ) from the underlying perturbative coupling a(Q 2 ) at high \Q 2 \ is proportional to a "sufficiently high" power of 
l/Q 2 . In Ref. [23], an analytic model was constructed for which Ai(Q 2 ) - a(Q 2 ) ~ (A 2 /Q 2 ) 3 at high |Q 2 | > A 2 . This 
was achieved by constructing the model-defining discontinuity function pi(a) = Im Ai(— a — ie) in such a way that 
it is equal to the discontinuity function pf t \a) = Im a(— a — ie) of the underlying pQCD coupling a(Q 2 ) at a > M 2 
(where Mq ~ 1 GeV), and in the low-a regime (0 < a < Mq) the unknown behavior of was parametrized 

by a single delta function. By adjusting the free parameters of the model, the aforementioned suppressed deviation 
Ai(Q 2 ) — a(Q 2 ) ~ (A 2 /Q 2 ) 3 at Q 2 > A 2 was achieved. 3 In Ref. [24j . this idea was continued, by employing in the 
low-energy regime (0 < a < Mq) the parametrization in terms of two delta functions for the discontinuity function 
pi(cr). In this way, the strongly suppressed deviation A\(Q 2 ) — a(Q 2 ) ~ (A 2 /Q 2 ) 5 at Q 2 > A 2 was achieved (while 
at the same time the correct value of r T was reproduced). Therefore, in this model the leading- twist contributions to 
physical quantities do not give power-suppressed terms of dimension D < 10 of ultraviolet origin. Hence, the model 
can be applied with OPE, without contradicting the ITEP interpretation of OPE, up to D — 8 terms. 

We wish to point out that our approach eliminates unphysical singularities from the QCD coupling while still 
preserving the salient features of the usual pQCD approach, among them the applicability of OPE (in the ITEP 
sense) and the related universality of the QCD coupling. On the other hand, the approaches of Refs. @, H, [26l - [28j 
follow a different line: by considering (various) specific timelike observables, they eliminate the unphysical singularities 
directly from the corresponding spacelike observables. This leads either directly 0, H, [HI or indirectly [23, HH to 
analytic QCD couplings whose nonunversality is reflected in observable-dependent modifications in the low-energy 
(low-cr) regime. 

In this work, we will apply the mentioned analytic model of Ref. [24| to the analysis of the r-decay data, in the V+A 
channel, using the OPE approach (with up to D = 6 terms) with Borel sum rules, along the lines of Refs. [29|, [3(| 
where the analysis was performed within the perturbative QCD in MS scheme. The program of the work is the 



following. In Sec. [II] we summarize the previously mentioned QCD analytic model 24j. In Sec. IIIII we recapitulate 
a powerful Pade-related resummation method |31| which is very natural and convenient to apply in evaluations of 
physical quantities within analytic QCD models. In Sec. IIVI we apply this resummation, within our model, in the 
evaluation of the (leading-dimension part of the) Adler function, i.e., the logarithmic derivative of the polarization 
operator n(Q 2 ) in the massless limit. We combine this evaluation, and the OPE expansion of the Adler function, with 
the Borel sum rules along rays in the complex plane of the Borel scales M 2 (as in Refs. [2!| [3(J; cf. also. Ref. [I2|). 
For the experimental input we use ALEPH results of 1998 for oj cxp (a) (oc Im n(— a — ie)) obtained from r-decay 
invariant- mass spectra for < a < m 2 [13] • We also discuss the obtained results and compare them with the results 
obtained when no resummation is performed in analytic QCD, and with the results in perturbative QCD in the 
renormalization scheme that is underlying the analytic QCD model (with and without resummation) and in MS 
renormalization scheme. Section |V] contains a summary of the obtained results and conclusions. 



II. 2-DELTA ANALYTIC QCD 

Here we recapitulate the analytic QCD model of Ref. [13], which contains a parametrization with two delta functions 
in the low- a regime of the discontinuity function pi(u) = Im Ai(— a — ie). Any analytic QCD model is defined, 
in principle, via its analytic spacelike coupling Ai(Q 2 ), which is the analytic analog of the perturbative coupling 
a(Q 2 ) = a s (Q 2 )/ir. Any other coupling A n (Q 2 ) (analog of the power a(Q 2 ) n ) can then be constructed from it, via 
the formalism introduced in Refs. [33j when n is integer and in Ref. [34| when n is general noninteger. 4 On the other 
hand, the coupling Ai(Q 2 ), analytic in the complex plane outside of the negative semiaxis Q 2 g C\(— oo,0], can be 



3 This condition was also achieved in the analytic model of Ref. |2T| . 

4 A formalism for construction of A n (Q 2 ) for general noninteger n, applicable only to MA, was developed and applied in Refs. [35l |. 
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expressed via the discontinuity function pi(a) = Im Ax(— a — ie) which is defined only for a > 0. The expression is 
a dispersion relation which follows from the application of the Cauchy theorem, the analyticity of Ai(Q 2 ), and the 
asymptotic freedom at large \Q 2 \ 

MQ 2 ) = -[ +C ° Vt?M- « 

7T Jo (<J + Q Z ) 

Due to the success of perturbative QCD at high energies, it is reasonable to assume that at higher a (a > Mq, 
where Mq > 1 GeV) the discontinuity function pi (a) coincides with the underlying perturbative function (a) 
1 r : i a{—a — ie). In the low-cr regime (0 < a < Mq) the behavior is unknown in its details, and is here parametrized 
with two delta functions 5 at lower values a = M 2 and Mf (0 < M 2 < Mf < Mq) 

2 

(%* V; 02) = tt]T f 2 A 2 S(a - Mf) + 6(<7 - M 2 ) x p^\a; c 2 ) (2) 

3=1 
2 

= ^E^ ,5 ( s -^) + e ( s ~ s «) xr i pt) ( s ; c 2) ■ ( 3 ) 

i=i 

In Eq. ([3]) dimensionless parameters were introduced: s — cr/A 2 , s 3 = M 2 /A 2 (j — 0,1,2), and r < f t \s;c2) — 

p^ t \a;c2) — Im a{Q 2 = —a — ie;c 2 ), c 2 being a scheme parameter. The Lambert scale A 2 (~ 10 _1 GeV 2 ) will be 
defined below. The underlying perturbative coupling a{Q 2 ) is taken, for convenience, in such renormalization schemes 
where it can be expressed explicitly as a solution in terms of the Lambert function W{z) 

a{Q 2 ; c 2 ) T 2 \ . (4) 

Here, Q 2 = \Q 2 \ exp(i</>); W-\ and W+i are the branches of the Lambert function for < <f> < +ir and — tt < <p < 0, 
respectively, 6 and the variable z involves the mentioned Lambert scale A 



1 f\Q 2 ^~ f3o/ci 



cie 



A 2 



exp(-i/3 o 0/ci) . (5) 



The explicit expression @ was presented in Ref. [40l |. It is the solution of the following (perturbative) renormalization 
group equation (RGE): 

da(Q 2 ]c 2 ) _ 2 [l + ( Cl - (c 2 /ci))a] 

ainQ 2 - - P ° a [1 - (c 2 / Cl )a] ' (6) 

Here, ^ = (1/4)(11 - 2n//3) and c x = fixlfo = (1/4)(102 - 38ri//3)/(ll - 2n//3) are universal constants. On the 
other hand, c 2 = /3 2 / 1 (3q is the free three-loop renormalization scheme parameter. The expansion of the above beta 
function (3(a) = da/d\nQ 2 gives 

/3(a) = -A)« 2 (l + cia + c 2 a 2 + c 3 a 3 + ...), (7) 

with the higher renormalization scheme parameters Cj = l3j//3 Q (j > 3) fixed by the value of c 2 : Cj = c^ _1 /cj~ 2 
(J>3)- 7 



5 We note that the function 

tt 7 M 2 (ct + Q 2 ) tt J M ? hr (<r + Q 2 ) 

is a Stieltjes function (and JV^hr ~ 10 _1 GeV is a QCD threshold scale). Approximating the discontinuity function p(<r) of any Stieltjes 
function f(Q 2 ) as a sum of delta functions is well motivated, cf. Refs. |36|. |37| . because it leads to approximating the Stieltjes function 
by near-to-diagonal Pade approximants. The latter must converge to the Stieltjes function when the order (i.e., the number of deltas) 
increases [3Sj | . An idea similar to Eq. but with one delta, was applied in Refs. [27| [28j directly to spectral functions of the vector 
current correlators. 

6 The functions W±i(z) are implemented in MATHEMATICA |«| by the commands ProductLog[±, z] . 

7 Further use of Lambert functions in QCD running is discussed also in Refs. [4ll444ll . 
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TABLE I: Values of the parameters of the considered two-delta anQCD model, under the restriction for the values of the 
pQCD-onset scale Mo = ^/soA.^■ 1 GeV < Mo < 1.5 GeV, and .4i(0) < 1.0. The corresponding Lambert scales A are for the 

central value of the QCD coupling parameter ai^(Mz) = 0.1184. 



C2 = /32//3o 


so 


Sl 


fi 


S2 


fi 


A [GeV] 


Mo 


-4i(0) 


-5.73 


25.01 


18.220 


0.3091 


0.7082 


0.6312 


0.231 


1.15 


1.00 


-4.76 


23.06 


16.837 


0.2713 


0.8077 


0.5409 


0.260 


1.25 


0.776 


-2.10 


17.09 


12.523 


0.1815 


0.7796 


0.3462 


0.363 


1.50 


0.544 



Application of the dispersion relation JT]) to the discontinuity function ([3]) gives the analytic coupling A\{Q 2 ) of 
the model 

where u — Q 2 /A 2 . The free parameters of the model (c2, A, sq, s\, f 2 , S2, /|) are fixed by various requirements. 
The scale A is fixed by requiring that the underlying perturbative coupling a(Q 2 ) reproduce the central value of the 
world average at Q 2 = Af§ as given in Ref. 0: a {MS) {M 2 z ) = 0.1184/tt. More specifically, the coupling a(Q 2 ;c 2 ), 
with a chosen value of c 2 and in the considered low-momentum regime of interest (Q < 2m c where n/ = 3), is first 
transformed to the exact four-loop MS scheme, then RGE-evolved up to Q 2 = M§ using the four-loop MS beta 
function and at the quark thresholds Q 2 hl = (2m q ) 2 using the three-loop matching conditions |46|. The Lambert 

scale A (at rif = 3) is then fixed such that the described procedure leads to the value a^ MS ^(M|) = 0.1184/7T. 
The four parameters Sj and f 2 (j = 1,2), at chosen values of so and C2, are determined as functions of sq by the 
earlier mentioned (ITEP-OPE motivated) requirement Ai(Q 2 ) — a(Q 2 ) ~ (A 2 /Q 2 ) 5 for \Q 2 \ > A 2 (these are four 
requirements, in fact). Subsequently, the value of the parameter sq (= Mq/A 2 ) is determined, for a chosen value 
of C2, by the requirement that the model reproduce the aforementioned experimental value of the V+A strangeless 
scmihadronic r decay ratio r T 0.203. Finally, the values of the parameter c 2 are varied in such a way that the 
resulting pQCD-onset scale M (~ 1 GeV) varies within a specific range. For phenomenological reasons, the preferred 
values of M should be below and not too close to the value of the r lepton mass (m T = 1.777 GeV). On the other 
hand, if Mq < 1.15 GeV, the values of A\{Q 2 ) become too large (Ai (0) > 1) and the model loses stability in the 
infrared. We choose the representative values of the scheme parameter ci such that the following conditions are 
fulfilled: 1 GeV < M < 1.5 GeV and Al(0) < 1- For the central value M Q = 1.25 GeV, we have the values of 
parameters given in the second line of Table |TJ The two border choices (Mo = 1 GeV; and Ai(0) = 1 with Mo ~ 1.15 
GeV) are also given in Table |U 8 

However, the world average of the coupling parameter [45[ has some uncertainty: a( MS )(M§ ) = (0.1184±0.0007)/7r, 
corresponding to a {ME \m 2 T ) nf=z = (0.3183 ±0.0057)/tt, i.e., A n/=3 = 0.336 ±0.010 GeV where A„ /=3 is the usual MS 
scale (at n/ = 3). This implies that the Lambert scale varies, A sa 0.260±0.008 GeV [a(m 2 , c 2 ) = (0.2905 ±0.0043)/?r], 
while all the (central) values of the dimensionless parameters of the model (02 = —4.76, so = 23.06, etc.) are 
unchanged. The predicted value of the V+A r-decay ratio then varies, r T = 0.203 ± 0.006, which is still compatible 
with the experimental value 0.203 ± 0.004 [H 

Having the analytic analog (a(Q 2 )) an = ./^(Q 2 ) in 2-delta QCD analytic model in Eq. ((5]), the analytization of 
higher integer powers (a™) an = A n is performed according to the construction in Refs. [33[ which is applicable to any 
analytic QCD model. We briefly present it below. The basic idea is to introduce the logarithmic derivatives 

5 »«<« 2 > S ^W< (-.= 1.2,.-.). (9) 

We note that a n+1 {Q 2 ) = a(Q 2 ) n+1 + 0(a n+2 ) by RGE da(Q 2 )/d\nQ 2 = j3(a(Q 2 )), where beta function (3(a) has 
the pQCD expansion as given in Eq. ([7]). Due to the linearity of analytization, it follows from (a(Q 2 )) an = Ai(Q 2 ) 
the relation ((9a(Q 2 )/(91n<3 2 )an = dA\{Q 2 ) / d\nQ 2 , and thus in general 

(S«+i(Q 2 )) an = A n+ i(Q 2 ) , (10) 



In Ref. l24ll . we included Mo = 1.00 GeV. However, in this case A.i(0) = 2.29, indicating instability in the infrared. 
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where 



A n+1 (Q Z ) = 



ffinl d(lnQ 2 ) 



2\n 



(n=l,2,...) 



(11) 



and where Ai(Q 2 ) is given in our case in Eq. ((SJ. An interesting aspect is that in virtually any analytic QCD 
model, including the present one, we have a clear hierarchy |^li(Q 2 )| > |^4 2 (Q 2 )| > 1-^3 (Q 2 ) I > ••■ not just for large 
\Q 2 \ > A 2 , but for any Q 2 , cf. curves in Fig. Ufa) (which are for the presented model and at Q 2 > 0). This suggests 



(b) 













(a) 


y, \ ', \ 








\ MS 






. a 










4,4,^ 












16 A, 














16 A, . ~ 











--16A3 



0.5 1.0 



5.0 10.0 



0.5 l.( 



5.0 10.0 



C 2 (GeV 2 ) 2 2 (OeV 2 ) 

FIG. 1: (a) The couplings Ai(Q 2 ), 4 x Az(Q 2 ) and 4 2 x Az{Q 2 ) of 2-delta analytic QCD model with the central input values 
(02 = —4.76; so = 23.06), as a function of positive Q 2 : 0.1 GeV 2 < Q 2 < 10.0 GeV 2 ; the rescaling factors 4 and 4 2 are for 
better visibility; included are the analogous pQCD couplings a, 4 x 0,2 and 4 2 x 013 in the same (Lambert) scheme, and the MS 

coupling a MS . (b) The couplings Ai{Q 2 ), 4 x A2{Q 2 ) and 4 2 x Az(Q 2 ), where the couplings A2 and ^3 are constructed by 
Eqs. (|22|I -(I23 | ) with truncation at (including) A4 term; included are the (Lambert) pQCD analogs a, 4a 2 and 4 2 a 3 . 

the following approach to the evaluation of any dimension-zero (D = 0) contribution T>(Q 2 ) of a massless spacelike 
observable, such as Adler function, in analytic QCD. Let the perturbation series (pt) of this quantity be 



2?(Q 2 ) P t = a(reQ 2 )^d„(re) a (reQ 2 )" +1 , 



(12) 



where /i 2 = kQ 2 is a renormalization scale, re ~ 1 being a fixed chosen dimensionless renormalization scale param- 
eter. 9 Before evaluating T>(Q 2 ) in analytic QCD, we reorganize the above series into the corresponding "modified" 
perturbation series (mpt) in logarithmic derivatives a n +i, Eq. §§§ 

00 

T> (Q 2 ) mpt = a(reQ 2 ) + £ d„(re) a n+1 {nQ 2 ) . (13) 
n=l 

This leads, after applying the analytization (|I0l) term-by-term, to the ("modified") analytic series (man) 

00 

2?(Q 2 )man - Al{nQ 2 ) + ^ d n {n) A n+ l{nQ 2 ) , (14) 

n=l 

which is the basic expression for evaluation of V(Q 2 ) in analytic QCD. 

Since the "mpt" series (fLU) is just a reorganization of the re-independent "pt" series (fT!Z1) . the "mpt" series V (Q 2 ) mp t 
is also re-independent. This then immediately implies, in conjunction with the recurrence relation <9a n (reQ 2 )/<91n re = 
— fiorian+i (kQ 2 ) [this being a direct consequence of the definition , the following set of differential relations between 
<2„(re): 



din. 



-dn (re) = nfiodn-i (re) (n=l,2, ...), 



(15) 



For k dependence of d n (K) coefficients, see the next Section IlIII Eqs. 1 1381 - d 101 ) . 
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where <io( K ) = do(n) = 1 by definition. Using these relations, it is straightforward to verify that the ("modified") 
analytic series T>(Q 2 ) man of Eq. (ITU) is K- independent. 10 

The coefficients d n (n) are obtained from <ifc(K)'s (k < n) in the following way. We first relate powers a n+1 and the 
logarithmic derivatives a n +i, at a given scale Q 2 (or: /i 2 = kQ 2 ) using the RGE relations in pQCD 11 



a 2 = a 2 + cia 3 + c 2 a 4 + • • • , (16) 
5 
2< 



5 

a,3 = a 3 + -cia 4 + • • • , 04 = a 4 + • • • , etc. , (17) 



and we invert them 



a 2 = a 2 - cia 3 + Qc 2 - c 2 ^ a 4 H , (18) 

5 

a 3 = a3 — -cia4 + ■ ■ • , a 4 = 04 + • • • , etc. (19) 

Replacing the relations (fT8|) - (fT9|) into the perturbation expansion (fl"2|) for T>(Q 2 ) we can read off the tilde coefficients 
d n (K) of the reorganized ("modified") expansions (TT^)) - (fm) 

= di(n) , d 2 (K) = d 2 {n) - c\d\(n) , (20) 
d 3 (n) = d 3 (n) - ^cid 2 (n) + Qc 2 - c 2 ^j di(n) , etc. (21) 

Applying analytization, Eqs. (|10p - (fTTj) , in relations (fT5 |) - (fT5|) term-by-term, we finally obtain the analytic analogs of 
integer powers, A n = (a") a „ 12 

A 2 = (a 2 ) an =A 2 - C!A 3 + Qc 2 -c 2 \m + --- , (22) 
A 3 = (a 3 ) an = A 3 - - Cl A4 + --- , Aa = (a 4 ) an = Aa + --- , etc. (23) 



This means that the "modified" analytic series (|14j) can be rewritten in the more usual form, in close analogy with 
the original perturbation series (fT2|) 

oo 

2?(Q 2 )an - AxinQ 2 ) + J2 <*n(«) A n+ i(nQ 2 ) . (24) 

71=1 

This series, being a reorganization of the K-independent series T>(Q 2 ) man of Eq. (Tl4|) . is therefore also n- independent. 
Several couplings A n in the present analytic QCD model, for low positive Q 2 , are presented in Fig.[ljb). 

In practice, in the expansions of T>(Q 2 ), Eqs. (TT2l and (IT3")) . we know exactly only a few first coefficients d n (n), up 
to (and including) n — n max , i.e., up to a™ max+1 . This implies that in the relations (fT6"|) - (fl~9"| and (|22p - (|2"3")l it is natural 
to perform truncations at (and including) the term ~ a" m " +1 (~ a nmax +i ~ -^•n max +i ~ -^n m ax+i)- For example, in 
the case of Adler function, n max = 3, the perturbation series and the reorganized series are truncated at a 4 (~ S4) 

3 

V(Q 2 -k) [ *1 = a( K Q 2 ) + Y,dn(n)a( K Q 2 ) n+1 , (25) 

71=1 

3 

2?(Q 2 ; k)W t = a(nQ 2 ) + £ d„( K ) a n+1 ( K Q 2 ) , (26) 



For k dependence of d n (n) coefficients, obtained upon integrating the relations iJTjjJl, see Eq. (I4.'it in the next Section lllll 
The first RGE is da(Q 2 )/d\nQ 2 = /3(a(Q 2 )), where beta function /3(a) has pQCD expansion given in Eq. iJTJ, with cj = fij/fio. The 
other RGE equations are obtained by applying further derivatives d/d In Q 2 to the first RGE. 

Expressions of An in terms of Ah's (k > n), which can be used via inversion to obtain An in terms of Ak's, for integer n and k, were 
given in the context of MA model of in Refs. U3, 133. |47| . 
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and the two analytic series (fl4t and (|24|) also become truncated, at A4 (~ A4) 

3 

^(Q 2 ; «)SL = -4i(«Q 2 ) + £ A l+ i(«Q 2 ) , (27) 

71=1 

3 

V(Q 2 \k)^1 = A 1 (KQ 2 ) + Y / d n ( K )A n+1 ( K Q 2 ) . (28) 

n=l 

If the relations (jTHJ) - ((TSJ) and (|2"2"l) - (|23l) are then truncated naturally, i.e., at (and including) ~ a 4 (~ 0,4 ~ A4 ~ Ai), 
it is straightforward to check that the truncated series (125|l is identical with (]26f) . and (f27f with (fl4j) . 

Due to the truncation, the above series are renormalization scale (/c) dependent. However, since the form of all the 
RGE relations of pQCD is maintained, by construction, also in analytic QCD under the correspondence a n H> A n 
and a n (-> A n , the truncated analytic series (|2"T|) and P5| have weak renormalization scale dependence of one order 
higher than the last included term 

— __ A(^)-^Q) aiT^ — ' ( 29 ) 

and this dependence is getting weaker when the number N = n max + 1 of terms in the truncated series increases 

dV(Q 2 ] /c)LJn r f „ n 2\ a t „rfi\ dV{Q 2 ; k)U , qn . 
am k In k 

The last relation ~ on the right side of Eq. (f3U| is valid as long as the couplings A n , appearing in V(Q 2 ; k)\^\ arc 
constructed via the linear combinations Eqs. (I22t - (|23|) with so many terms that the last term Am included there has 
M > N. The derivative of the "man" truncated series on the left-hand side of Eq. (f5U|) is in fact only one term 



dV(Q 2 ;n 



[N] 

man 



9 In n 



-PoN<1n-i{k)An+i ( k Q 2 ) , (31) 



as can be explicitly checked by using the recursive relation 13 dA n +i{nQ 2 ) /dhx k — —Po(n + l)A n +2(nQ 2 ) and the 
differential relations (p~5)) between the coefficients d n (n) which are a consequence of k independence of the full "man" 
series The derivative of the "an" truncated series on the right-hand side of Eq. (15tJ|) is in general a finite 

linear combination of terms Ar+\(kQ 2 ) with R> N, the number of the terms of this combination depending on the 
level of truncation made in the construction of A n 's in the relations (|22 |) -(f23 |) . One of the benefits of using analytic 
QCD is that this residual unphysical dependence [~ An+i(kQ 2 ) ~ An+i(kQ 2 )] is getting weaker with increasing N 
irrespective of the physical momentum scale Q 2 (in stark contrast with truncated series in pQCD), because of the 
aforementioned hierarchy |„4i(Q 2 )| > |^2(Q 2 )| > |^3(Q 2 )| > • ■ ■ which is valid at any Q 2 (not just: \Q 2 \ > A 2 ). The 
described method of analytization of integer powers a n was constructed in Refs. [33], and was extended in Ref. [13] 
to the case of terms with noninteger powers a v and of terms of the form a v ln fc a. 

If the analytization of higher powers a n were performed in the naive nonlinear way, (a ra ) an = A", the renormalization 
scale dependence of the truncated series would in practice not decrease with the inclusion of more terms in the series, 
but would in general even increase, because in such a case the derivatives d/d In k of the truncated series include 
complicated nonperturbative contributions such as 1/(Q 2 ) K , cf. Appendix C of Ref. [22J (second entry). The full 
(naive) analytic power series T>(Q 2 ) aps = Ai(kQ 2 ) + J2^Li d n {K){Ai{nQ 2 )) n+1 is not K-independent, basically because 
the powers Ai(nQ 2 ) n do not fulfill RGE relations analogous to those of a(nQ 2 ) n . 

The RGE of the coupling A\{Q 2 ) has now, by construction, formally the same form as the RGE of pQCD coupling 
a(Q 2 ) where we replace a™ \-> A n 



dA^Q 2 
d\nQ 2 



(= -p A 2 {Q 2 )) = -PoA 2 {Q 2 ) - PiA 3 (Q 2 ) - foA^Q 2 ) (32) 



In Sec. II VI we will present the curves of the resulting Ai(Q 2 ) and of the underlying Lambert pQCD coupling a(Q 2 ) 
on the <5 2 -contours in the complex plane which are needed for the evaluation of the Borel sum rules there. 



13 This being a direct consequence of the definition illl . 
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In the extraction of the experimental value of the strangeless V+A decay ratio r T (AS — 0, m q — 0) CX p = 0.203 ± 
0.004, [3 [l9j], the contributions of the higher dimension (D — 2,4,6,8) chirality-violating terms (i.e., nonzero 
quark mass effects) were subtracted. These latter terms were estimated to be Sr T (AS — 0,to u ^ ^ 0) = (— 5.8 ± 
1.4) x 10~ 3 (cf. also Appendix B of Ref. [22| ). The chirality-nonviolating contributions were not subtracted from 
the mentioned r T (AS = 0,m q = 0) cxp = 0.203 ± 0.004. Among the chirality-nonviolating D > 2 contributions, the 
only possibly nonnegligible one (cf. Ref. i30j) is the D = 4 contribution from the gluon condensate, (Sr T )( a< 3G) = 
(ll/4)a.g(m^.)(aGG)/m^.. However, in our evaluation of r T (AS — 0,m q = 0) in the analytic QCD we assumed that 
the gluon condensate contribution is negligible, i.e., r T (AS = 0,m q = 0) was evaluated as the D = contribution 
only, and was required to achieve the value 0.203. Later in this article, we will deduce that the gluon condensate in 
analytic QCD has similar values as in pQCD, i.e., (aGG) « 0.005 GeV 4 , which then gives us the contribution to r T 

(64 V+A A a l^!l 2 ( m 2) J_ ( aGG) w L4 x 10 -4 (33) 
V T J (aGG) 4 my 

where we replaced a 2 (m 2 )/7r 2 = a(m 2 ) 2 (« a 2 (m 2 )) by .A 2 (ra 2 ) w 0.01, cf. Fig. Ufa). Thus we can conclude a 
posteriori that our neglecting of the higher dimension chirality-nonviolating terms in the OPE of r T was justified in 
our evaluation of r>, where the latter evaluation contributed significantly to the fixing of parameters of the (2-delta) 
analytic QCD model. 



III. PADE-RELATED RESUMMATION OF ADLER FUNCTION IN ANALYTIC QCD 



In this Section we summarize an evaluation method for massless spacelike physical quantities, which we apply 
to the evaluation of the dimension D = contribution of Adler function, T)(Q 2 ;D = 0) = T>{Q 2 ). This method 
was developed some time ago for pQCD evaluations [H, EH and is a generalization of the diagonal Pade (dPA) 
resummation method in pQCD We will first apply dPA to the mentioned (D = 0) Adler function T>(Q 2 ). 

The perturbation series of this quantity is known up to the fourth term, Eq. (|25[) . where /i 2 = kQ 2 is the (squared) 
spacelike renormalization scale (k ~ 1), and the truncated series has a residual /x 2 -dependence due to truncation. 
Coefficients d n (l) = dj(l,MS) (n = 1,2,3), in MS scheme and at the renormalization scale /x 2 = Q 2 , were obtained 
in Refs. 51-53], respectively 



= + 0.691772 /3 n , 
= -27.849 + 8.22612 /3 + 3.10345/3, 



Mi) 

d 3 (l) = 32. 727 - 115. 199 f3 + 49. 5237 (3 2 + 2.18004 



o • 



(34) 

(35) 
(36) 



The light-by-light contributions were not included in these coefficients, since they are zero when the number of effective 
quark flavors is rif — 3 (then: j3o = 9/4). The value rif = 3 is used in the evaluation of T>(Q 2 ) because the relevant 
energies in the analysis of the next Section are \Q 2 \ < m 2 (m 2 w 3.2 GeV 2 < (2m c ) 2 ~ 6.5 GeV 2 ). 

We are interested in evaluation of Adler function in 2-delta analytic QCD model of Sec. |TT]and in the corresponding 
pQCD with the same renormalization scheme. Therefore, we have to transform first this expansion from the MS 
renormalization scheme to the new (Lambert) scheme of 2-delta analytic QCD model: C2 = — 4.76;c3 = c^/ci; etc., 
cf. second line of Table HI and Eqs. (HJ)-© which define the corresponding pQCD and the renormalization scheme. 
The scheme invariance of the perturbation expansion (|12p then implies that the coefficients in the new (Lambert) 
scheme are 



di(l) = di(l) 



d 2 (l)=d 2 (l)-( C2 -c 2 ) 



ds(l) = d 3 (l)-2d 1 (l)(c 2 -c 2 )- -(c 3 -c 3 ) 



(37) 



where the bars denote the values in MS scheme. The new expansion coefficients c? J (/i 2 /<5 2 ), in Lambert scheme and 
at the renormalization scale /i 2 = kQ 2 , are then 

di(/s) = di(l)+A)hi(/«) , 

2 \ qfcl„fc 



da(«) = d 2 (l)+^l ~ ]^ln"(K)da_ fc (l)+^iln(K) , 
fc=i ^ ' 

d 3 (n) = d 3 (l)+Y / (l)l3^n( K )d^ k (l)+(3 1 
fe=i ^ ' 



21n(/s)di(l) + -A) ln'(«) 



/3 2 ln(«) 



(38) 
(39) 

(40) 
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where these relations were obtained from the renormalization scale invariance of the perturbation expansion (| 12[) . The 
resulting truncated perturbation expansion V(Q 2 ; k)^, Eq. (|25|) . is then used as the basis for the evaluation of the 
D = Adler function T>(Q 2 ) in 2-delta analytic QCD model. Due to truncation, it has (unphysical) renormalization 
scale dependence, dT>(Q 2 ; K)^/dhiK ~ a 5 , and the same is true for the corresponding analytic truncated series 

V{Q 2 -n) [ X Eqs. tf28j), i.e., 9X>(QV)S/91n/t~ A 5 , cf. Eq. (29]). 
The dPA-resummed result can be written in two equivalent ways 14 



[2/2]p(a( M 2 )) - 



1 + Fix + F 2 x 2 



(41) 

x— a(/i 2 ) 
X 



ai- h a 2 



1 + U\X 1 + 



(42) 

x — a(/i 2 ) 



where aj. + as = 1. In Ref. |5G | it was shown that this approximant is independent of the renormalization scale 
fi 2 — kQ 2 (i.e., independent of k) used in the original truncated series (|2"5|) if the RGE-running is at the one-loop 
level. Building on this idea, in Refs. [48| this approach was extended so that the ^-independence of the resummed 
result was exact. 15 For this, the truncated perturbation series T>(Q 2 ; n) p l , Eq. (|2"5)l . in powers of a(/i 2 ), was first 

reorganized into the truncated series 2?(Q 2 )[^ pt , in logarithmic derivatives <x n +i defined in Eq. ©. The factor in front 
of this definition was chosen so that a\ = a and a n+ i = a n+1 + 0(a n+2 ) for n > 1. We recall that at one-loop level 
a n +i = a n+1 (in general: a n+ i ^ a n+1 ). 

The reorganized truncated (modified) perturbation series T>(Q 2 ; K)|^ pt is given in Eq. (|26[) . where the new coefficients 
dj(K) are related to the original coefficients dj(K) in Eqs. (|20l) - (|2"Tj) . and Cj = f3j/fto are the coefficients of the beta 
function, Eq. ©. The coefficients dj(n) can be regarded as "the one-loop parts" of the original coefficients dj(n), 
because they have the simple one- loop type of renormalization scale dependence (involving only the /?o coefficient). 
Namely, upon integrating directly the differential relations (|15[) . which were obtained on the basis of n independence 
of the full series T)(Q 2 ) mp t of Eq. (|T3l) . we obtain the explicit form of k dependence of d n (K.) 

d n (K) = d n (l) + E ( ) Po ln fc (^)d„- fc (l) , (43) 

k=l ^ ' 

where we recall that k is the dimensionless renormalization scale parameter (k = [i 2 /Q 2 ), and da — do = 1. The 
procedure for the construction of the generalization of dPA method consists now in the following. In the modified 
truncated series T>(Q 2 ; n)^j pt , Eq. (f26)) . we replace, in the one- loop sense, the logarithmic derivatives by the powers 
a n (/j, 2 ) i — y au(fj, 2 ) n , and obtain a truncated power series of a new quantity T>(Q 2 ) 

3 

V(Q 2 ; k)^ = a u ( K Q 2 ) + £ d n {n) a u ( K Q 2 ) n+1 , (44) 

n=l 

where o,u(kQ 2 ) — au(ji 2 ) is the coupling RGE-evolved from a(Q 2 ) to a(ji 2 ) by the one-loop beta function (—/3 a 2 ) 

au{^Q 2 ) = (n2 , ■ (45) 

The full (untruncated) series V{Q 2 ) = au(nQ 2 ) + Y^ = i d n (n) au{nQ 2 ) n+1 is k- independent, as can be easily checked 
by using the differential relations fp~5|) . Now we apply to the truncated power series of T>(Q 2 ) pt the [2/2] Padc 



We recall that in the case of a four- term power series l|25|l . the diagonal Pade is [2/2] (a), which is by definition the ratio of two quadratic 
polynomials in a = a(fj, 2 ) such that T> — [2/2] (a) ~ a 5 . 

Since the physical quantity X>(Q 2 ) is /^-independent, this extended rcsummation, having the same property, is expected to approximate 
better the (unknown) full expression for T>(Q 2 ). 
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approximant, as was performed in Eqs. (|41|) - (j42]) in the case of T>(Q 2 

x(l + Exx) 



[2/2] 5 {au{nQ 2 )) = 



F x x- 



(46) 



~ X „ X 

oli - — h a 2 - 



(47) 

x=a\i (kQ 2 ) 



1 + Ui(k)x 1 + U 2 (n)x 

= 5i aui^Q 2 ) + a 2 a u (n 2 Q 2 ) . (48) 
Here, 5i + a 2 = 1- Going from Eq. (|4"Tf to (|4"8")l . we used the relation (l4"5"j) and denoted the new scales as 

Q 2 = KjQ 2 . = k exp , (j = l,2). (49) 

The resulting approximant for the original truncated power series (I25[) is then obtained by simply replacing the 
one-loop pQCD coupling an in Eq. (|4"8"|) by the full pQCD coupling a (= a s /ir) 

G l v 2] (Q 2 ) = "I a(^iQ 2 ) + a 2 a(n 2 Q 2 ) , (50) 

This method of construction can be applied in a completely analogous way when the number N of known perturbation 
terms in the series of T>(Q 2 ) pt is any even number N — 2M (N = 2,4, 6, . . .), 16 leading to the approximant 

M 

gjf/^Q^f^a^Q 2 ), (51) 

where . . + 5m = 1. In Rcf. [48| it was proven that the result is exactly independent of the original renormalization 
scale n 2 = kQ 2 . In the proof in Ref. [48| it was demonstrated that each weight coefficient a.j and each scale coefficient 
Kj is separately independent of the renormalization scale parameter k; for a somewhat less formal and more intuitive 
proof, see Appendix [XI here. The K-independent coefficients a.j and Kj = Q 2 /Q 2 are also Q 2 -independent since they 
are dimensionless. In Ref. [48| it was also proven that the approximant fulfills the basic approximation requirement 
required of any resummation approximant 

V{Q 2 ) - <4 M/M] (Q 2 ) = 0(a 2M+1 ) = 0(a 2M+1 ) . (52) 

The approximant (|50p . although theoretically attractive, turned out not to work well within pQCD. The reason for 
this is that one of the two scales in Eq. (|50|) . e.g. Qf, is lower than Q 2 and often brings the pQCD coupling a(Q\) 
close to the Landau singularities, making thus the evaluation of low-momentum physical quantities T>(Q 2 ) in pQCD 
unreliable. For example, in MS scheme, Adler function (with nf = 3) gives K\ <C 1 (cf. also footnote [25] in Sec. II VI) . 

However, in Ref. [3l| this method was revived and applied in analytic QCD frameworks, where no such problems 
of Landau singularities appear. It was demonstrated in Ref. [3l| that the approximant (|50[) should be applied with 
the same weights and the same scales as in pQCD also in analytic QCD where the analytic truncated power series of 
the physical quantity has the form (|28|) analogous to Eq. (|25|) and the modified truncated analytic series (|27|) has the 
form analogous to Eq. (|26|) . 17 The resummed result is completely analogous to Eq. ([50)) 

^ /2] (Q 2 ;an) = Si A^Q 2 ) + a 2 Ai{n 2 Q 2 ) , (53) 

with dij and Kj obtained by construction in Eqs. (|4"?| and (l4"9l and thus /i 2 - and Q 2 -independent. The applicability 
of the approximant (f53|) is based on the fact, proven in Ref. [3lj, that it also fulfills (the analytic analog of) the basic 
approximation requirement, i.e., 

V{Q 2 ) - 4 2/2l (Q 2 ;an.) ~ MQ 2 ) ~ A 5 (Q 2 ) . (54) 
For the analogously constructed q\^^ m ^ (Q 2 ; an.), the above difference becomes ~ A 2 m+i (~ A 2 m+i), cf. Ref. [3l| . 



When N = 2, the method reduces to the (two-loop) effective charge method of Refs. [54l . 

The couplings Aj+i(n 2 ) are obtained from Ai(fJ?) in complete analogy with Eq. J9|, according to Eq. The construction of 

^4 n (/x 2 )'s, as a linear combination of these quantities, is given in Eqs. | |22I) - ||23| |. and was explained in detail in Refs. [33, 34]. 
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We recall that in pQCD we regard V = V pt = T> mpt , and in analytic QCD V = T> an = T> man , where these series 
quantities are written in Eqs. (fT2"j) - (fll?)) and in Eqs. ([23)1 and (TT4")) , respectively. 

In the renormalization scheme of 2-delta analytic QCD model (Lambert scheme central choice: C2 = —4.76, Cj = 
ci^ 1 /c[~ 2 for j > 3), where the three coefficients dj(n) at general renormalization scale parameter k (= fi 2 /Q 2 ) are 
obtained via Eqs. (|3"4"l) -(f40" |) . the described formalism gives us the following values of the scale coefficients Kj and 
weights OLf 

ki (= Ql/Q 2 ^ = 0.1689 , Si = 0.6418 , 

k 2 (= Qj/Q 2 ^ = 3.1656 , 5 2 = 0.3582 . (55) 

These quantities are each exactly independent of the choice of the original renormalization scale parameter k (= fi 2 /Q 2 ) 
in the original expansion coefficients (|38l) - (|40|) . as proven in Ref. [48( and in Appendix [A] and as can be checked also 
numerically by starting with the construction of these quantities from c?j (k')'s at a different k' . 

The massless Adler function T>(Q 2 ) = T>(Q 2 ; D = 0), which is a logarithmic derivative of the leading-dimension part 
of the polarization operator (correlator) of hadronic currents, will play a central role in the next Section in the analysis 
of the Borel sum rules involving invariant-mass spectra of the r lepton decay, applied within 2-delta analytic QCD 
model described in the previous Section. We will thus apply the method of resummation described here, Eq. (|53[) . for 
the evaluation of the Adler function at complex Q 2 (\Q 2 \ ~ m 2 ). 



IV. ANALYSIS OF TAU DECAY DATA WITH BOREL SUM RULES 



The idea of sum rules in r decay physics could be summarized as an application of the identity 

/ daf(-a)u} e * p (a) = -m <t dQ 2 f(Q 2 )U th {Q 2 ) , (56) 

JO J|Q2| =CTo 

where the contour integration on the right-hand side is in the counterclockwise direction, f(Q 2 ) is an analytic function 
in the Q 2 complex plane, and w(cr) is the spectral function of the polarization function Tl(Q 2 ) of hadronic currents 

w(cr) = 2vr Im II(Q 2 = -a - it) . (57) 

The identity (|56|) is obtained by applying the Cauchy theorem to the function /(Q 2 )n(Q 2 ) and taking into account the 
analytic properties of the physical polarization function Ii{Q 2 ) as required by the general principles of quantum field 
theories. We recall that in pQCD-evaluated (pQCD+OPE) polarization function n(Q 2 )th these analyticity properties 
are in general not respected, because of the Landau singularities of the pQCD coupling a(Q 2 ). This means that in 
pQCD, when we replace on the left-hand side of Eq. (|56|) oj cxp (a) H> cjth(c), the identity in general ceases being valid. 
In analytic QCD no such conceptual problems appear, the theoretically evaluated n(Q 2 ) t h automatically respects the 
analyticity properties on which the sum rule (|56p is based. 

In this work we are interested in the nonstrange V+A channel of r decays. As a consequence, the polarization 
function is a sum of functions 

n(Q 2 )^4 1) (Q 2 ) + n«(Q 2 ) + ng ) (Q 2 ) . (58) 

These functions appear in the polarization operators IT^, (q) which are correlators of the (nonstrange) charged hadronic 
currents 



KM = i J d^xe^iTJ^J^) 



= (q»q» - 9»vq 2 )W{Q 2 ) + q^Uy (Q 2 ) , (59) 
where : J = V, A; J M = uj^d ( J = V) , J^= vrf^d (J = A) . (60) 

For more details on these points, we refer to Refs. [29|, [3(|. On the left-hand side of the rum rule (|56|) the experimental 
spectral function w OX p(c) is used, obtained from the measured r-decay invariant-mass spectra for < a < m 2 . 18 On 



The data published by ALEPH Collaboration are in Refs. [l7l , ll8ll . The experimental bands represented by the left-hand side of Eq. (156 b . 
for f(Q 2 ) = exp(Q 2 /M 2 )/M 2 , are taken here from Figs. 4 and 5(a), (b) of Ref. [23], which in turn were obtained from the ALEPH 
Collaboration data of 1998, Ref. [TJ]. 
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the right-hand side of Eq. (|56l) . the theoretically evaluated polarization function H(Q 2 ) appears, which can be evaluated 
with the OPE 

n(Q 2 )th = ~ ln(Q 2 /^ 2 ) + n(Q 2 ; D = 0) + £ |M (l + C n a(Q 2 )) • (61) 

The D — 2 operator term (i.e., n = 1 term) comes from the nonzero values of the current masses of u and d quarks, 
it is negligible and is neglected here. For the evaluation of the contour integral on the right-hand side of Eq. ([56]). it 
is convenient for us to perform integration by parts, resulting in 

jH d*f(-a)u, exp (a) =-± jT_^ ^V Adl (Q 2 ) [HQ 2 ) - H-so)] | Q2=soexpW) , (62) 

where 2?Adi(Q 2 ) is full Adler function, i.e., 

2? Adl (Q 2 ) = -2^^^ (63) 



d\nQ 

and the function T is any function satisfying the relation 

dT{Q 2 ^ 



l+P(Q 2 ;^ = 0) + 2^^I^(l+C n a(Q 2 )) , (64) 



dQ 2 

The dimension zero {D = 0, or leading-twist) terms are related by 



f(Q ) ■ (65) 



(V(Q 2 ) =) V(Q 2 ;D = 0) = -2^ ^±^ . (66) 



The perturbation expansion of the D — massless and strangeless Adler function T>(Q ) = T>(Q ' , D = 0), cf. Eq. (|25p . 
is now known up to ~ a 4 order with the coefficients dj(l) (j = 1,2,3) in MS scheme and at renormalization scale 
/i 2 = Q 2 given in Eqs. ([34"1) - (|36|) for the here relevant momentum regime \Q 2 \ < m 2 (i.e., for nj — 3). We transformed 
the expansion from MS renormalization scheme to the new (Lambert) scheme of 2-delta analytic QCD model (c 2 = 
—4.76; C3 = cfi/ci; etc., cf. second line of TableHJ according to the relations (f37|). and the resulting coefHcients dj{n), 
at an arbitrary renormalization scale parameter k = /i 2 /Q 2 , are given in Eqs. (|38[) - (|40[) . The resulting truncated 
analytic series, in 2-delta analytic QCD, is given in Eqs. (|27 |) - (|28| . The latter two expressions are identical because the 
construction of the analytic analogs (a") an = A n was performed as a truncated linear combination of the logarithmic 
derivatives ^4fe's of A\ [cf. Eq. pip] according to Eqs. (|2"2")) - ([2U)) , with the last included term in those linear combinations 
being Aa- On the other hand, the resummed expression was obtained in Eq. (|53")) in 2-delta analytic QCD and Eq. ([50)1 
in pQCD, with the RGE-invariant values of the scale and weight coefficients Kj and 5y given in Eqs. (|55|) . We refer 
for more details to Sections ITll and Hill 

The basic idea of the Borel sum rules is to choose, in the sum rule relation (|56|) [or: (|62])]. for the function / an 
exponential function [2^, [3(| 

f(Q 2 ) = ^ exp(Q 2 /Af 2 ) , F(Q 2 ) = exp(Q 2 /M 2 ) , (67) 

where M 2 are, in principle, arbitrary complex scales. Other choices of function /, in the context of sum rules, have 
also been used in the literature, cf. Refs. [31, [EB-lHl] ■ The integrals in the sum rules ([56]), ([62]) . with the choice ([67]) . 
become Borel transforms 19 B(M 2 ), and the sum rule acquires the form 

B cxp (M 2 ) = B th (M 2 ) , (68) 



We choose in the sum rules H56II and 1 1621) for the upper integration bound the largest possible value, i.e., uq = m^. If co is significantly 
below , it has been argued in the literature that in such a case the duality- violating effects become important and have to be accounted 
for, cf. Ref. 
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where 



B cxp (M 2 ) = 
B th (M 2 ) = B(M 2 ;D = 0) + 2tt 2 Y^ 



exp(-o-/M 2 )w cxp (cr) , 
(O 



2,i 



n>2 



(n - 1)! (M 2 )« ' 



and the D = part is the following contour integral: 

B(A/ 2 ;D = 0) = — [ d(f>V(Q 2 =m 2 T e l ' i ';D = 0) 

27T 



exp 



M 2 



exp 



IT 2 



(69) 
(70) 

(71) 



The Borel transform suppresses the D > 4 terms by (n — 1)! factor (n = D/2) 20 and suppresses the high energy tail 
of Wcxp(c) where the experimental errors are larger. For the low-energy regime \Q 2 \ < 1 GeV 2 , it does not provide 



suppression. In Refs. [29|, |3fJ it was argued that the D = part of the Borel sum rule can be reliably calculated 
within pQCD only for \M\ 2 > 0.8-1 GeV 2 , due to the (unphysical) Landau singularities of the pQCD coupling a(Q 2 ). 
In analytic QCD we do not have this problem. 

Eq. (jTTj) implies that we have to evaluate the Adler function D(Q 2 ; D — 0) along the contour Q 2 = m 2 exp{i(j>). In 
particular, if resummation is performed, Eqs. (I53[) and (|50[) imply that the relevant coupling parameters in the Adler 
function are Ai(Kjm 2 exp(i(j))) and a(Kjm 2 exp(i^)), with k\ and «2 as given in Eq. (1551) . The real and imaginary parts 
of these couplings are presented in Figs. [2] (a), (b). We can see from these Figures that the analytic coupling A\ and the 
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FIG. 2: (a) Real and (b) imaginary part of the analytic coupling Ai(Q^) and of the underlying pQCD coupling a(Qj) (j = 1, 2) 
for the contour Q 2 = m 2 . exp(j0). At the higher scales (Q 2 ), the two couplings are indistinguishable. The scheme parameter 
value is C2 = —4.76 (Lambert scheme, cf. the second line of TableU}. 



underlying pQCD coupling a differ from each other appreciably only at the low scales, i.e., at Q\ = 0.1689 m 2 exp(i^). 

On the other hand, if not performing the resummation in analytic QCD, the truncated series (|27j) and thus the 
logarithmic derivatives Aj+\(ii 2 ) have to be evaluated; and in the underlying pQCD, the powers of a([i 2 ) have to be 
evaluated, cf. Eq. (|2"5j) . In all such cases, the otherwise arbitrary renormalization scale fi 2 = kQ 2 will be chosen with 
k = 1, i.e., Q 2 = m 2 exp(i</>). The real and imaginary parts of the couplings A\(Q 2 ) and a(Q 2 ), as well as of the MS 
coupling a(Q 2 ; MS), as functions of the contour angle <f>, are presented in Figs. [3] (a), (b). We can see that the analytic 
coupling Ai(Q 2 ) and the underlying pQCD coupling a(Q 2 ) are almost indistinguishable at such scales \Q 2 \ — m 2 . 
In fact, also the corresponding logarithmic derivatives cij+i(Q 2 ) and Aj+i(Q 2 ) are very close to each other. 21 This 
means that the two "modified" truncated approaches (|2"o| and (|27p. with /i 2 = Q 2 (= m 2 exp(i<f>)), give us results very 
close to each other. While in analytic QCD we will truncate the series in terms of logarithmic derivatives, Eq. (|27l) . 



We checked that the terms C„a(Q 2 )(0 2 n)/{Q 2 ) n in the OPE expansion Il61|l give negligible contributions and were thus ingored in the 
Borel sum rules J68l - (f70T . 

Since we have in 2-delta analytic QCD the relation Ai(Q 2 ) — a(Q 2 ) ~ (A 2 /Q 2 ) 5 (for \Q 2 \ > A 2 ), repeated application of d/d In Q 2 to 
this relation gives Aj+i(Q 2 ) - a J+ i(Q 2 ) ~ (A 2 /Q 2 ) 5 . 



14 




FIG. 3: (a) Real and (b) imaginary part of the analytic coupling Ai(Q 2 ) and of the underlying pQCD coupling a(Q 2 ) for the 
contour Q 2 = m 2 exp(i<j!>), with C2 = —4.76. The two couplings are practically indistiguishable. Included is also MS coupling 
(c 2 w 4.47, c 3 « 20.99) along the contour. 



we will, however, in pQCD apply the truncation to the power series, Eq. (|2"5j) (with /x 2 = Q 2 ), instead. Therefore, 
due to this somewhat different kind of truncation, the difference will be appreciable, but small nonetheless, between 
the truncated analytic and the truncated pQCD approaches. 

In order to separate or isolate terms of various dimensions in the Borel sum rule (|S5)) - (jT0)) . the Borel transform is 
evaluated along fixed chosen rays in the complex M 2 plane [2!| [3(| 

M 2 = |Af| 2 exp(iV) , 0.68 GeV 2 < \M\ 2 < 1.5 GeV 2 . (72) 

For example, when if> = ir/6 (ip = tt/4), the real part of the Borel transform contains no D — 6 (D = 4) term 22 
because Re(e l7T / 2 ) = 



R Cj B cxp (|M| 2 e"/ 6 ) = RcB(|M| 2 e -/ 6 ;I9 = 0) + 7r 2 ^ , (73) 



RcB cxp (\M\ 2 e™/ 4 ) = R c B(\M\ 2 e™/ 4 -,D = 0)-Tr 2 ^2^- , (74) 

V2|M |° 



The D = 4 and D — 6 operators can be expressed in terms of condensates [2C 



1 

1287T 2 



(Of = ^(aG^G^) , (75) 



/n( V + A )\ ~ L/,0n I- \2 tna\ 
{°6 ) ~ — ^~ a ^> ' ( 76 ) 

where a = a s /i:. The approximation (1761) for (Og V ' + ' 4 '') is obtained after factorization (vacuum saturation) assumption 
of various 4-quark condensate contributions, and is expected to be valid with not better than 20-30% accuracy [20j,|29|. 

The results of our analysis for the ip = ir/6 are given in Fig. 0] Using 2-delta analytic QCD model described in 
Sec. HI1 and the resummation method of Sec. lIIII for the evaluation of the D = part, Eq. (I7T|) . we obtain for the gluon 
condensate the values 

(aGp„Gp V ) = (0.0055 ± 0.0040 (cxp) ± 0.0025 (oth) ) GeV 4 = 0.0055 ± 0.0047 GeV 4 . (77) 

For the parameters of the analytic QCD model we used the central values of the parameters as determined in Ref. [24j |. 
i.e., the second line of Table U here. The central condensate value 0.0055 GeV 4 was obtained by adjusting the gluon 
condensate value in such a way that the theoretical curve (the solid line in Fig. [4]) gives the minimal deviation 



We include in the OPE sum H61M64H only terms up to n = 3, i.e., dimension D = 2n = 6. 
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0.8 1.0 1.2 1.4 

|M| 2 /GeV 2 

FIG. 4: Real part of the Borel transform, Eqs. l|68 ]l -l(70 )) and 1(73)1. for M 2 = \M\ 2 exp(#) with V = tt/6. The band with dark 
vertical stripes represents the experimental data (see footnote 1 18[) . The solid line is the best theoretical curve of the resummed 
expression for D — Adler function T>(Q 2 ) in 2-delta analytic QCD model (anQCD resummed), obtained by the standard 
minimization with respect to the central experimental values. It corresponds to the central value of the gluon condensate 
(aGG) — 0.0055 GeV 4 . The light grey band represents the variation of this curve when the standard minimization is applied 
with respect to the values of the upper and lower bounds of the experimental band, and gives the variation of the gluon 
condensate 0.0055 ± 0.0040 GeV 4 . 



from the central experimental values, 23 in the standard minimization procedure, using 40 equidistant points over the 
depicted |M| 2 -interval. The resulting deviation parameter was x 2 = 6.40. An interesting feature of the resulting 
best theoretical curve (solid line) is that it remains within the experimental band in almost the entire considered 
interval of |M| 2 , only slightly surpassing the upper bound of the experimental band at the lowest values of |M| 2 
(|M| 2 w 0.68 GeV 2 ). 

The uncertainty S(aGG) cxp = ±0.0040 GeV 4 was obtained by applying the same minimization procedure to the 
upper bounds and lower bounds of the experimental band, respectively. The resulting two theoretical curves then 
define the upper and the lower border of the light grey band in Fig. |4j This band is thus the prediction of the method 
(resummed analytic QCD), without the other uncertainties included. 

The other uncertainty S(aGG) ot h = ±0.0028 GeV 4 was obtained as the variation of the resummed analytic QCD 

prediction when the QCD coupling parameter value is varied in its world average interval [H], ai MS) (M 2 ) = 0.1184=F 
0.0007; and when the scheme parameter c-i of the analytic model is varied according to Table HI C2 = — 4.76 +2 'g7. 
At the end of Sec. HH these two variations in the model are discussed in more detail. 24 The calculations give us for 
uncertainties of the value of the gluon condensate coming from these two effects: 8(aGG) as =±0.0025 GeV4 > and 
S(aGG) C2 =|o;qoo4 GeV 4 . Adding in quadrature then gives us the "other" uncertainty 6{aGG) ot h ~ ±0.0025 GeV 4 . 

If we adjust the gluon condensate values, by the same standard minimization method (40 points) with respect to 
the experimental values, in various evaluation approaches [anQCD resummed; anQCD (truncated); Lambert pQCD 
resummed; Lambert pQCD (truncated); MS (truncated)], we obtain the following values for the gluon condensates, 



At a given value of |M| 2 , the central experimental value was considered to be the arithmetic average of the upper and the lower bound 
value of the experimental band at that |M| 2 . 

ai MS) (M|) = 0.1184 corresponds to: a< MS ) (m 2 )„, =3 = 0.3183/vr (and to the standard MS scale: A n/=3 = 0.336 GeV); and in Lambert 
scheme with C2 = —4.76 to a(m 2 ) = 0.2905/7T (and to Lambert scale A = 0.260 GeV). On the other hand, varying the renormalization 
scale fi 2 = kQ 2 = Km 2 exp(i</>) (<-> varying re) in the coefficients dj(re) of the expansion 1 1251) does not change at all the resummed results 
for the D = Adler function V(Q 2 ) appearing in the contour integral 1711 . as argued in Sec. I II 1 1 
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FIG. 5: Real part of the Borel transform, Eqs. (55) -® and lffi|l. for A/ 2 = |M| 2 exp(it/>) with V = tt/4. The band with 
dark vertical stripes represents the experimental data (see footnote 1 18[). The best theoretical curve (solid line) is obtained by 
standard minimization of the resummed analytic QCD results with respect to the central experimental values. The light grey 
band is obtained by applying the same minimization to the upper and lower bounds of the experimental band. The dashed 
curve is the prediction (of the resummed analytic QCD) when (Og +A ^) = 0. See the text for details. 



and the corresponding deviation parameter \ 2 '- 



(0.0055 ±0.0047) 


GeV 4 


x 2 


= 6.40 , 


(anQCD resummed) , 


(0.0104 ±0.0058) 


GeV 4 


x 2 


= 9.60 , 


(anQCD truncated) , 


(0.0056 ± 0.0049) 


GeV 4 


x 2 


= 7.42 , 


(Lambert resummed) , 


(0.0122 ±0.0046) 


GeV 4 


x 2 


= 9.83 , 


(Lambert truncated) , 


(0.0059 ±0.0049) 


GeV 4 


x 2 


= 7.43 , 


(MS truncated) . 



(78) 



The uncertainties given above were obtained by adding in quadrature the experimental uncertainty <5(aGG) C xp = 
±0.0040 GeV 4 and the uncertainty from <5ai MS) (Af§) = ^0.0007 [=> 5a [ ™ S) (m 2 T ) nf=3 = ^0.0057] which is 
S(aGG) as = ±0.0025, ±0.0025, ±0.0028, ±0.0023, ±0.0028 GeV 4 , for the respective five methods above. In the 
two analytic QCD approaches (anQCD resummed, anQCD truncated) the aforementioned uncertainty coming from 
c 2 = -4.761q is S(aGG) C2 =+o;qqo4 GeV 4 and +0.0012 GeV 4 , respectively, and was also added in quadrature. We 
checked that the gluon condensate values change by less than 2% if the number of points in the standard minimiza- 
tion is decreased from 40 to 28; the \ 2 values decrease, but the values of ratios of \ 2 between different methods are 
maintained. 

It is interesting that the evaluation with the pQCD approach in Lambert scheme (resummed) and in MS scheme 
(not resummed), 25 give similar results for the gluon condensate as in anQCD resummed (and with only somewhat 
higher \ 2 )- 

However, the latter is not the case for the curves when ip = 7r/4 [cf. Eq. ([73)1 ] . The results for the case ip = ir/4 
are presented in Fig. [5j The central theoretical curve (solid line) represents the resummed analytic QCD result 



We did not apply the resummation approach of Sec. IIIII to D = Adler function in pQCD in MS scheme. Namely, it turns out that 
MS is such a scheme which requires for the Adler function that the lower of the two invariant scales be very low: |£Ji| <C 1 GeV (while 
\Q\ = m T pb 1.78 GeV), leading to severe problems with the Landau singularities. In the Lambert scheme, on the other hand, we obtain 
|Qi| si 0.73 GeV and \Q 2 \ w 3.16 GeV. 
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when the condensate value (0^ /+A ' > ) is obtained with the standard minimization (with 40 equidistant points), in 

complete analogy with the ip — tt/6 case of Fig. 01 This central value is (0 < ^ /+A " > ) = —0.5 x 10~ 3 GcV 6 . The 
resulting deviation parameter is \ 2 = 10.79. The best theoretical curve (solid line) only slightly surpasses the 
upper bound of the experimental band at low \M\ 2 < 0.7 GeV 2 , remains within the experimental band in the 
interval 0.7 GeV 2 < \M\ 2 < 1.0 GeV 2 , and is situated slightly below the lower bound of the experimental band for 
1.0 GeV 2 < |M| 2 . It is interesting that the theoretical curve with (0 { / +A) ) = (dashed line) is, in comparison, not 
singificantly worse: it is situated slightly below the lower bound of the experimental band for 0.9 GeV 2 < \M\ 2 , and 
its deviation parameter is x 2 — 12.65. 

The uncertainty S(0^ +A ^) cxp = ±0.9 x 10~ 3 GeV 6 is obtained by applying the same minimization procedure to 
the upper bounds and lower bounds of the experimental band, respectively. The resulting two theoretical curves then 
define the the light grey band in Fig. [5J in analogy with the light grey band of Fig. 0J 

The other uncertainty of D = 6 condensate is estimated as coming from the coupling parameter a s and scheme 
parameter c 2 uncertainties, which are: <5(0^ +j4) ) Qs = ±0.4 x 10" 3 GeV 6 and 5(0 { ^ +A) ) C2 = ±0.4 x 10~ 3 GeV 6 . 
Adding in quadrature this gives us 5{0^ +A " l ) o t\ l — ±0.6 x 10~ 3 GeV 6 . We thus obtain the following estimate for 
D = 6 condensate of the V+A channel: 

(0^ +A) ) = (-0.5 ± 0.9 (oxp) ± 0.6 (oth) ) x 10~ 3 GeV 6 « (-0.5 ± 1.1) x 10~ 3 GeV 6 . (79) 

The standard minimization gave for the central value a negative number close to zero (-0.5 x 10~ 3 GeV 6 ). The result 
(|T9"]) corresponds approximately to the factorized quark condensate values of Eq. (fT6f 



a <??> 2 * 7^<°6 > « (-3- ± 7.) x 10- 5 GeV 6 . (80) 



81 l (V+A) 

128tt 2 K 

We note that the estimates (j79 )) -([80 l) are not incompatible with nonnegative values. Further, they are not incompatible 
even with the following positive values of a(qq} 2 extracted from the V-A sum rules of Ref. [30j 26 : 



(Oj^ A) ) = -(4.4 ±0.6) x 10" 3 GeV 6 (81) 



a{qqf 



64tt 2 x 1.33 



(O, 



(V-A), 



(4.7 ±0.6) x 10" 5 GeV 6 . (82) 



It is interesting that the methods other then the resummed analytic QCD approach give us for the central value of 
D = 6 condensate significantly more negative values. Specifically, if we adjust in these methods the D — 6 condensate 
value by the standard minimization (with 40 equidistant points) to the experimental values, we obtain the following 
values of the condensate, and of the \ 2 fitting parameter: 

(0^ V+A) ) = (-0.5 ± 1.1) x 10~ 3 GeV 6 , X 2 = 10.79 , (anQCD resummed) , 

= (-1.9 ± 1.3) x 10~ 3 GeV 6 , \ 2 = 15-75 , (anQCD truncated) , 

= (-1.8 ± 0.9) x 10~ 3 GeV 6 , \ 2 = 12-23 , (Lambert resummed) , 

= (-2.3 ± 0.9) x 10~ 3 GeV 6 , x 2 = 16-15 , (Lambert truncated) , 

= (-1.8 ± 0.9) x 10~ 3 GeV 6 , X 2 = 12-38 , (MS truncated) . (83) 

The uncertainties above include, in quadrature, the experimental uncertainty <5(Oe)exp = ±0.9 x 10~ 3 GeV 6 , and the 
uncertainty from <5ai MS) (M 2 ) = ±0.0007, which is S(0 6 ) as = ±0.4 x 10~ 3 GeV 6 for the two analytic QCD methods, 
and ±0.3 x 10~ 3 GeV 6 for the three pQCD methods. In the two analytic QCD approaches (anQCD resummed, anQCD 
truncated) the uncertainty coming from c 2 = -4.76t^;^ is 6(0 6 )c 2 = (toj) x 10~ 3 GeV 6 and x 10" 3 GeV 6 , 

respectively, and was also included in quadrature. The results of all these methods, with the exception of our central 
method (anQCD resummed ), a re incompatible with nonnegative values of {0^ +A ^) . This agrees also with the analyses 
of r-decay data in Refs. [Ill [57]], where finite energy sum rules were applied (in pQCD±OPE approach) in MS scheme. 
The results (|83p thus suggest that the factorization assumption leading to the relation (f76")l fails to predict correctly 
even the sign of the condensate, i.e., that it fails much more severely than by 20-30% mentioned in Ref. [2{|. On the 



V-A sum rules are more adequate to determine the quark condensate a{qq) 2 , because in such a case the factorization approximation, 
in contrast to the V+A case Eq. H76H . does not involve subtractions of large terms. 
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other hand, the properly resummed analytic QCD (+ OPE) approach gives us the result (l79l) and (|80|) . suggesting 
that the relation ([75)1 does not necessarily fail so severely. In fact, the resummed analytic QCD gives for the choice 

(Og 1/+j4 ' ) ) = the value \ 2 = 12.65, which is quite similar to the values y 2 = 12.38, 12.23 obtained by the minimization 
in MS pQCD and Lambert resummed pQCD, respectively (for the central value (Og V+ " 4 ' ) ) ~ —1.8 x 10~ 3 GeV 6 ). 
When we increased the number of (equidistant in |M| 2 ) points in the standard minimization procedure from 40 

to 60, the values of (0^ +A * > ) in Eqs. (|83|) changed very little, and did not affect the digits displayed there; and the 
corresponding values of x 2 increased proportionally, numerically by about 49-50%. 

Finally, when -0 = 0, the Borel transform is real and includes, in principle, both condensate contributions (D = 4, 6). 
It is presented in Fig.[S] The result of the analytic QCD (with resummed D = Adlcr function), for the central values 
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FIG. 6: The Borel transform for real M 2 = \M\ 2 (ip = 0). The band with dark vertical stripes represents the experimental 
data (see footnote I18|) . The light grey band represents the variation of the resummed 2-delta analytic QCD model prediction 
(anQCD resummed) when D = 4 and D — 6 condensates vary in quadrature around their central values 0.0055 GeV 4 and 
-0.5 x 10 -3 GeV 6 , Eqs. £77]) and JTJJ. The central solid curve is for the central values of both condensates. The dashed curve 

is for the choice (0^ +A ^) = (and (aGG) = 0.0055 GeV 4 ). See the text for details. The results of other methods are given 
as points, for their corresponding central values of D — 4 and D = 6 condensates. 



of the condensates, Eqs. (|77[) and (|79l) . is represented by the central solid curve which is within the experimental band 
in the entire interval of M 2 . The effects of variation of the values of the two condensates, Eqs. (|TT[) and (f79|) . are also 
presented in Fig. El by the light grey band; the total variations of the two condensates, S(02n) (n — 2, 3), are added 
in quadrature [cf. Eq. ([70)) ]: 



B{M ) light grey 



= B{M 2 ;D = Q) + 2it 2 ^ 



(0 2 n) 



n=2,3 



1)! (M 2 



±2tt z 



S(Q 4 )\ 2 (5(0 6 )\ 2 
M 4 J V 2!M 6 J 



1/2 



(84) 



The dashed curve is the result of the choice (0 { ^ +A) ) = (and (aGG) = 0.0055 GeV 4 ). The results of other methods, 
with their corresponding central values of the condensates, Eqs. (|78p and (|53")) . are included as points in Fig.[Bl These 
results are close to the lower bound of the experimental band for all M 2 . The resulting \ 2 values (with 40 points) , with 
respect to the central experimental results (with -0 = 0) are: \ 2 = 4-44 (anQCD resummed); x 2 — 13.44 (anQCD); 
X 2 = 7.34 (Lambert resummed); x 2 — 13.84 (Lambert); x 2 = 7.58 (MS). It is interesting that the dashed line (anQCD 
resummed with (0 { ^ +A) ) = and (aGG) = 0.0055 GeV 4 ) gives x 2 — 7.02, which is even slightly less than MS (with 
its own best central values: (0 { / +A) ) = -1.8 x 10~ 3 GeV 6 and (aGG) = 0.0059 GeV 4 ). 
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V. CONCLUSIONS 

In this work we reanalyzed the Borel sum rules along the rays in the M 2 -complex plane, for the strangeless 
V+A r-decay invariant-mass spectra of ALEPH Collaboration of 1998. The analysis was performed here with the 
analytic QCD model developed in Ref. (24|. In this model, the analytic coupling A\(Q 2 ) has no unphysical (Landau) 
singularities, it predicts the correct measured value of V+A strangeless decay ratio of r lepton, and at high squared 
momenta it differs from the underlying perturbative coupling a(Q 2 ) (= a s (Q 2 )/ir) by terms ~ (A 2 /Q 2 ) 5 , allowing us 
to keep the ITEP school interpretation of the OPE expansion up to (and including) dimension-eight (D = 8) terms. 
The analysis involves the evaluation of an integral of the D = Adler function T)(Q 2 ; D = 0) along the contour \Q 2 \ = 
m 2 . We evaluate this function by applying a generalization of the diagonal Pade resummation, the method being 
exactly invariant under variation of the renormalization scale and very convenient for application within analytic QCD 
framework (the latter because no Landau singularities appear). By comparing the experimental ALEPH data of 1998 
with the theoretical evaluations for the Borel transform along the argument ray M 2 = \M\ 2 exp(iijj) with ip = 7r/6, we 
extract with the standard minimization for the (D = 4) gluon condensate the values (aGG) = (0.0055±0.0047) GeV 4 , 
cf. Fig. 0] and Eq. (fTT)) . Furthermore, considering the ray ip — 7r/4, we analogously extract for D = 6 condensate 
the values (0 { / +A) ) = (-0.5 ± 1.1) x 1(T 3 GeV 6 , cf. Fig. [5] and Eqs. (JTHJ) - dHni) • This is not incompatible with the 
theoretical expectation (0^ +A ^) oc a(qq) 2 > based on the factorization (vacuum saturation) approximation. Using 

the obtained central values of the two condensates, (aGG) = 0.0055 GeV 4 and (Oq V+A ^) = —0.5 x 10~ 3 GeV 6 , the 
Borel transform on the real positive axis ("0 = 0) gives the theoretical curve which remains within the experimental 
band in the entire analyzed MMnterval (0.68 GeV 2 < M 2 < 1.5 GeV 2 ), cf. Fig.© 

We also compare these results of the resummed analytic QCD approach with those of the not resummed analytic 
QCD approach and with those of perturbative QCD (pQCD) approaches in two schemes: Lambert scheme of the 
aformentioned analytic QCD model (c-i = = —4.76, c 3 = c^/ci, etc.), and in MS scheme (c2 ~ 4.47, C3 w 20.99). 

It turns out that in the Lambert scheme pQCD resummed approach and in the MS pQCD nonresummed approach 
the extracted central values of the gluon condensate (aGG) are similar, 0.0056 and 0.0059 GeV 4 , respectively, and 
the x 2 values are only somewhat higher (by about 15%). On the other hand, in the not resummed pQCD approach 
in Lambert scheme, and in the not resummed analytic QCD approach, the extracted values of the gluon condensate 
are significantly different, with the central value of about 0.0122 and 0.0104 GeV 4 , respectively, and % 2 values are by 
about 50% higher than in the resummed analytic QCD approach. On the other side, the values of D — 6 condensate 
extracted (with i\> = 7r/4) in these approaches are neg ative, (0 { ^ +A) ) (-2± 1) x 10~ 3 GeV 6 , suggesting a complete 

failure of the vacuum saturation approximation for (0^ +A ^) (cf. also Refs. [56l[57j on this point). Furthermore, these 
approaches give us, for their corresponding central values of D = 4 and D = 6 condensates, at ip — the results which 
are situated close to the lower edge of the experimental band for all considered M 2 . 

Our result for the gluon condensate, (aGG) = (0.0055 ± 0.0047) GeV 4 , is similar to the results of Refs. l30l l56l . 
In Ref. [23], the value (0.006 ± 0.012) GeV 4 was obtained from the V+A channel of r-decay data, and in Ref. [3fJ 
the value (0.005 ± 0.004) GeV 4 from combination of the latter (with ALEPH 1998 data) and of the charmonium sum 
rules. In Ref. |56|, where ALEPH 2005 data were used, the approach of weighted finite energy sum rules with the 
contour-improved perturbation theory gave approximately the values (0.008 ± 0.005) GeV 4 for the gluon condensate 
when the standard MS QCD scale was 27 A n/= 3 = 0.30-0.35 GeV, the higher values when A = 0.3 GeV and the lower 
values when A = 0.35 GeV. The original work on the sum rules, Ref. [20I ]. obtains a higher estimate, of 0.012 GeV 4 
from charmonium physics. QCD-moment and QCD-exponential moment sum rules for heavy quarkonia give higher 
values, (aGG) (0.022 ± 0.004) GeV 4 and (0.024 ± 0.006) GeV 4 , respectively, Refs. [Hj]. On the other hand, in 
Ref. [19| a combined fit to the V+A r-decay data primarily extracts the value of a s and, as a byproduct, obtains the 
gluon condensate value (aGG) = (—0.015 ± 0.003) GeV 4 which differs significantly from the aforementioned values. 

The interpretation of our results can be summarized in the following. The nonexistence of the unphysical (Landau) 
singularities in the running coupling of the (2-delta-parametrized) analytic QCD framework of Ref. [24[ leads to 
a reasonable degree of consistency between the theory and the experiment in the Borel sum rule analysis of the 
strangeless V+A channel of r-decay physics. To achieve this, it is crucial to perform, in addition, an efficient Pade- 
related and renormalization scale independent resummation of the Adler function, a method which is very convenient 
for applications within analytic QCD frameworks precisely because of the absence of the Landau singularities. 



27 This corresponds to a{ MS \m^;n f = 3) between 0.2986 and 0.3259. 
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Appendix A: Renormalization scale independence of parameters in the resummed expression 

In this Appendix we show that the dimensionless parameters otj and Kj (j — 1,2), obtained via the construction in 
Eqs. (pH |) - ([4"9")) . are exactly independent of the (dimensionless) renormalization scale parameter k (= /j 2 /Q 2 ). 

We recall that the truncated series V(Q 2 ; k)^3 of Eq. (|44[) is made of d n (n) coefficients which are constructed 
from the coefficients (4(k) of the original power series (|25l) via the relations (|2T)|) - (l2"Tj) . and whose k dependence is 
a consequence of k independence of the full perturbation series V(Q 2 ) pt and V(Q 2 ) mpt of Adler function V(Q 2 ) (or 
any spacelike observable) [Eqs. (fT2 > ([T5 ]) ]. Further, the full series V(Q 2 ) = a u {nQ 2 ) + a u (KQ 2 ) n+1 is 

K-independent, because of the differential relations (1X5]) . 

Let us choose two different renormalization scale parameters, k and k' , and define the following notations 

x = au(nQ 2 ) , x' = au{n'Q 2 ) , (Al) 
u(k) = /3 ln/t . (A2) 

The running of the coupling an is one- loop, Eq. (|45j) . With the notation (|A2|) . the quantities Uj(k)'s appearing in 
the construction in Eq. (j47|) [see also Eq. (j49|) ] can be written as 



u j( K ) = u ( K j/ K ) ( = Po l n ( K j/ K )) ■ (A3) 
It is straightforward to verify from Eq. (|45p for the one-loop running coupling that x and x' are related by 

(A4) 



1 + u(k' j k)x ' 

and that the simple fractions appearing in Eq. (|47|) can be reexpressed as 



l + iij(n)x \ 1 + u^Kj j ' k)x J 1 + u(Kj/n')x' 



(A5) 



Further, it is straightforward to verify that the truncated series T>(Q 2 )^ of Eq. (|4"4")l , while being K-dependent due to 
truncation, at the two renormalization scale parameters k and k' differ from each other by ~ x 5 (~ x' 5 ) terms 

P K ,(x') - P K (x) ~ x' 5 (~ x 5 ) , (A6) 

where the two truncated series of Eq. (|44j) with n and k' are denoted as P K (x) and P K >(x') which are quartic polynomials 
in x and x' , respectively 

3 

P K (x) = V{Q 2 -n)f t ^x + Y J dn{n)x n+1 , (A7) 

n=l 
3 

P K ,{x') = V(Q 2 ; = x' + dnWW n+1 ■ (A8) 



Namely, the relation (|A6[) follows from the fact that both P K (x) and P K /(x') differ from the full ^-independent series 

V{Q 2 ) = a u (nQ 2 ) + £^=1 d n {n) a u ^Q 2 ) n+1 by terms ~ x 5 (~ x' 5 ). 

On the other hand, Eq. (|A5j) implies for the sum of simple fractions appearing in the construction of the method, 
Eq. ijlTjl. the following identity: 

2 x 2 x' 

^ ^l + UiWa; = ^ ^'TT^T^y ' (A9) 
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The left-hand side of Eq. (|A9I) is, by construction, the diagonal Pade [2/2]p K (x) of the polynomial P K (x) [= V(Q 2 ; k)^], 
cf. Eqs. (j46])-(|47]). 28 Therefore, Eq. l|A"9)l in conjunction with Eq. (TA"6|) implies that the right-hand side of Eq. (TA"9|) is 
the diagonal Pade [2/2] p k , (x') of the polynomial P K i{x') [= T>(Q 2 ; «/)pt ]. More explicitly, we have 

g^- TT^V "^^ = (g^T+^- p - (a!) ) +(^)-^(x'))- 5 -' 5 ■ (Aio) 

The difference in the first pair of parentheses on the right-hand side of Eq. (|A10[) is ~ x 5 because the sum of simple 
fractions there is, by construction, the diagonal Pade [2/2]p K (x) of the polynomial P K (x); the difference in the second 
pair of parentheses is ~ x 5 due to the relation (|A6|) . Eq. (IA10[) means that the sum of the simple fractions on the 
left-hand side there is the diagonal Pade [2/2}p k , (x') of the polynomial P K <(x'). 

However, since P K i {x 1 ) is a polynomial whose coefficients are entirely independent of k (they are only ^'-dependent), 
this leads to the conclusion that the coefficients of the Pade [2/2)p , (x 1 ) appearing on the right-hand side of Eq. (|A9[) 
are k- independent, i.e., that oij and Uj(nj/n') are ^-independent. This means that Kj and ctj, which were obtained 
via the construction in Eqs. (|44l) - (|49[) . are /{-independent. This concludes the demonstration. 

The above argument can be almost literally repeated for the general case of N = 2M terms (M = 1,2,3,.. .) in 
the original truncated series (|23|) and thus in the related truncated series T>{Q 2 ] k')^^ . The encountered diagonal 
Pade's are now [M/M], and the decomposition in simple fractions is a sum of M simple fractions. 
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